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Abstract 

The //z/SSM provides a solution to the //-problem of the MSSM and explains 
the origin of neutrino masses by simply using right-handed neutrino superfields. 
Given that R-parity is broken in this model, the gravitino is a natural candi- 
date for dark matter since its lifetime becomes much longer than the age of the 
Universe. We consider the implications of gravitino dark matter in the /ii^SSM, 
analyzing in particular the prospects for detecting gamma rays from decaying 
gravitinos. If the gravitino explains the whole dark matter component, a grav- 
itino mass larger than 20 GeV is disfavored by the isotropic diffuse photon back- 
ground measurements. On the other hand, a gravitino with a mass range between 
0.1 — 20 GeV gives rise to a signal that might be observed by the FERMI satel- 
lite. In this way important regions of the parameter space of the /ii^SSM can be 
checked. 



1 Introduction 



The "/i from u" Supersymmctric Standard Model (/iz/SSM) was proposed in the htera- 
ture [1-3] as an ahernative to the Minimal Supersymmetric Standard Model (MSSM). 
In particular, it provides a solution to the //-problem [4] of the MSSM and explains the 
origin of neutrino masses by simply using right-handed neutrino superfields. 

The superpotential of the /iz/SSM contains, in addition to the usual Yukawas for 
quarks and charged leptons, Yukawas for neutrinos L 0'^, terms of the type u^HdHu 
producing an effective /i term through right-handed sneutrino vacuum expectation 
values (VEVs), and also terms of the type v^v'^i)'^ avoiding the existence of a Goldstone 
boson and contributing to generate effective Majorana masses for neutrinos at the 
electroweak scale. Actually, the explicit breaking of R-parity in this model by the 
above terms produces the mixing of neutralinos with left- and right-handed neutrinos, 
and as a consequence a generalized matrix of the seesaw type that gives rise at tree 
level to three light eigenvalues corresponding to neutrino masses [1] . 

The breaking of R-parity can easily be understood if we realize that in the limit 
where Yukawas for neutrinos are vanishing, the v'^ are just ordinary singlet superfields, 
without any connection with neutrinos, and this model would coincide (although with 
three instead of one singlet) with the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) where R-parity is conserved. Once we switch on the neutrino Yukawa cou- 
plings, the fields v'^ become right-handed neutrino superfields, and, as a consequence, 
R-parity is broken. Indeed this breaking is small because, as mentioned above, we have 
an electroweak-scale seesaw, implying neutrino Yukawa couplings no larger than 10~^ 
(like the electron Yukawa). 

The latter also implies that processes violating lepton number that might wash-out 
any baryon asymmetry present in the model would be suppressed. Notice also that 
electroweak baryogcnesis could work in this model in a similar way to the case of the 
NMSSM [5]. Actually, the fact that in the //z/SSM there are three singlets instead of 
one like in the NMSSM, should in principle give more freedom to be able to obtain 
more easily electroweak baryogcnesis. The detail conditions for baryogcnesis in this 
model are presently under study [6]. 

Since R-parity is broken in the yUZ/SSM, one could worry about fast proton decay 
through the usual baryon and lepton number violating operators of the MSSM. Never- 
theless, the choice of -R-parity is ad hoc. There are other discrete symmetries, like e.g. 
baryon triality which only forbids the baryon violating operators [7]. Obviously, for all 
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these symmetries R-parity is violated. Besides, in string constructions the matter su- 
perfields can be located in different sectors of the compact space or have different extra 
U{1) charges, in such a way that some operators violating i?-parity can be forbidden 
[8], but others can be allowed. 

Several recent papers have studied different aspects of the /xz/SSM. In [2], the param- 
eter space of the model was analyzed in detail, studying the viable regions which avoid 
false minima and tachyons, as well as fulfill the Landau pole constraint. The structure 
of the mass matrices, and the associated particle spectrum was also computed, paying 
special attention to the mass of the lightest Higgs. In [9], neutrino masses and mixing 
angles were discussed, as well as the decays of the lightest neutralino to two body {W- 
lepton) final states. The correlations of the decay branching ratios with the neutrino 
mixing angles were studied as another possible test of the /xz/SSM at the LHC. The 
phenomenology of the //i/SSM was also studied in [10], particularized for one and two 
generations of right-handed sneutrinos, and taking into account all possible final states 
when studying the decays of the lightest neutrahno. Possible signatures that might 
allow to distinguish this model from other R-parity breaking models were discussed 
qualitatively in these two works [9, 10]. In [11], the analysis of the vacua of the //i/SSM 
carried out in [2] was completed, obtaining that spontaneous CP violation through 
complex Higgs and sneutrino VEVs is possible. Neutrino physics and the associated 
electroweak seesaw mechanism was also studied. It was shown how the experimental 
results can easily be reproduced and explained why the mixing patterns are so different 
in the quark and lepton sectors. All the results were discussed in the general case with 
phases. 

On the other hand, when it!-parity is broken, the lightest supersymmetric particle 
(LSP) is no longer stable. Thus neutrahnos [12] or sneutrinos [13], with very short fife- 
times, are no longer candidates for the dark matter of the Universe. Nevertheless, if the 
gravitino is the LSP its decay is suppressed both by the gravitational interaction and 
by the small R-parity violating coupling, and as a consequence its lifetime can be much 
longer than the age of the Universe [14]. Thus the gravitino can be in principle a dark 
matter candidate in R-parity breaking models. This possibility and its phenomenolog- 
ical consequences were studied mainly in the context of bilinear or trilinear R-parity 
violation scenarios in [14-22]. In [16, 18, 19, 21] the prospects for detecting gamma rays 
from decaying gravitinos in satellite experiments were also analyzed. In this work we 
want to discuss these issues, gravitino dark matter and its possible detection in the 
FERMI satelfite [23], in the context of the piuSSM. 
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2 The /iz/SSM 



The superpotential of the /^.z/SSM introduced in [1] is given by 

where we take Hj = {HlHa). Hi = (^+,^°), Qj = Lj = (/>,,e,), = 

1, 2, 3 are family indices, a,b — 1,2 are SU{2)l indices with €±2 — 1, and Y, A, k are 
dimensionless matrices, a vector, and a totally symmetric tensor, respectively. 

Working in the framework of supergravity, the Lagrangian >Csoft is given by: 



+ Cab 



+ (A,r,),,i/^L«z>; + c.c. 



+ 



-eay>{AyX), H^H'^ + ^^^^1 + c.c. 

^ (^Mg A3 A3 + M2 A2 A2 + Ml Ai Ai + c.c.) . 



(2) 



In addition to terms from £soft5 the tree-level scalar potential receives the D and F term 
contributions also computed in [1] . In the following wc will assume for simplicity that all 
parameters in the potential are real. Once the clcctrowcak symmetry is spontaneously 
broken, the neutral scalars develop in general the following VEVs: 



{i>i)^Vi, (yf)-^^- 



(3) 



For our computation below we are interested in the neutral fermion mass matrix. 
As explained in [1,2], neutralinos mix with the neutrinos and therefore in a basis where 
(5°, W^, Hd, Hu, I'Ri, J^Li), one obtains the following neutral fermion mass terms 
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The above 10 x 10 matrix, Eq. (5), is of the seesaw type giving rise to the neutrino 
masses which have to be very small. This is the case since the entries of the matrix M 
are much larger than the ones in the matrix m. Notice in this respect that the entries 
of M are of the order of the electroweak scale while the ones in m are of the order of 
the Dirac masses for the neutrinos [1, 2]. 

At low energy the free parameters of the //z/SSM in the neutral scalar sector are [2]: 
\i, Kijk, niHa, rnn^, mi,,, mvf., Ax^, A^^^^, and A,y^^. Strong upper bounds upon the 
intergenerational scalar mixing exist, so in the following we assume that such mixings 
are negligible, and therefore the sfermion soft mass matrices are diagonal in the flavour 
space. Thus using the eight minimization conditions for the neutral scalar potential, 
one can eliminate the soft masses mn^, itlh^, itt-z^ ^i^d m^rc in favour of the VEVs Vd, 
Vu, Vii and z/f. On the other hand, using the Standard Model Higgs VEV, v ~ 174 
GeV, tan/3, and z/j, one can determine the SUSY Higgs VEVs, Vd and f„, through 

= + + ^i- t^ViS consider as independent parameters the following set of 
variables: 



\-i i^ijki tan^, I'iji'i, Ax^, A^..^, A^.^ 



(8) 



It is worth remarking here that, because of the minimization conditions, the VEVs of 
the left-handed sneutrinos, z/j, are in general small, of the order of Dirac masses for the 
neutrinos [1]. Then, since Vi « Vd,Vu we can define the above value of tan /3 as usual, 
tan/3 = ^. 
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We will assume for simplicity that there is no intergenerational mixing in the pa- 
rameters of the model, and that in general they have the same values for the three 
families. In the case of neutrino parameters, following the discussion in [2,11], we 
need at least two generations with different VEVs and couplings in order to obtain the 
correct experimental pattern. We choose Y^-^ ^ Y^^ = ^1^3 S'^d I'l 1^2 = i^3- Thus the 
low-energy free parameters in our analysis are 

A, K, tan^, i/i, 1/3,1/"=, Ax, A^, A^ , (9) 

where we have defined X = \i, k, = Km, = z/f. Ax = Ax^, A« = A^^^^, A^, = Ay... 
Actually, we have checked that with Y^,,^ = Y^,^ ^ 2 Fj,^ ~ 10"^ and U2 = 1^3 ^ 
2 z/i ~ 10~^ GeV, the observed neutrino masses and mixing angles arc reproduced. 
As explained in detail in [11], this result is obtained so easily due to the peculiar 
characteristics of this seesaw, where R-parity is broken and the relevant scale is the 
electrowcak scale. 

The soft SUSY-breaking terms, namely gaugino masses, M12.3, scalar masses, 
uIq -c gc, and trilinear parameters, Au^ii,e, are also taken as free parameters and spec- 
ified at low scale. 



3 Gravitino dark matter 



Let us now show that the lifetime of the gravitino LSP is typically much longer than 
the age of the Universe in the //i/SSM, and therefore it can be in principle a candidate 
for dark matter. In the supergravity Lagrangian there is an interaction term between 
the gravitino, the field strength for the photon, and the photino. Since, as discussed 
above, due to the breaking of R-parity the photino and the left-handed neutrinos are 
mixed, the gravitino will be able to decay through the interaction term into a photon 
and a neutrino [14]^. Thus one obtains: 

r(«,/,^ET''.)-5^|c..Plf . (10) 

^ Other possible decay modes such as gravitino decay into a and a charged lepton, or into a 
Z'^ and a neutrino [19] are not relevant in our case, since we will obtain below that a gravitino mass 
smaller than 20 GeV is convenient in order to hilfiU experimental constraints. Neither we consider 
the possibility that the gravitino might in principle decay to singlet Higgs-neutrino if the Higgs is 
sufficiently light. 
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where m3/2 is the gravitino mass, Mp = 2.4 x 10^^ GeV is the reduced Planck mass, 
and |^/7;/p determines the photino content of the neutrino 

3 

lU^A"^ = ^\Niicosew + Ni2smew\'^ ■ (H) 

i=l 

Here Nn {Ni2) is the Bino (Wino) component of the i-neutrino. 
The Ufetime of the gravitino can then be written as 



,2^ -1 



If \Uj„\'^ ~ 10~^^ — 10~^^ in order to reproduce neutrino masses, as we will show below, 
the gravitino will be very long lived as expected (recall that the lifetime of the Universe 
is about 10^^ s). 

For the gravitino to be a good dark matter candidate we still need to check that 
it can be present in the right amount to explain the relic density inferred by WMAP, 
0,DMh^ ~ 0.1 [24]. With the introduction of inflation, the primordial gravitinos are 
diluted during the exponential expansion of the Universe. Nevertheless, after inflation, 
in the reheating process, the gravitinos are reproduced again from the relativistic par- 
ticles in the thermal bath. The yield of gravitinos from the scatterings is proportional 
to the reheating temperature, Tr, and estimated to be [25] 

where nig is the gluino mass. As is well known, adjusting the reheating temperature 
one can reproduce the correct relic density for each possible value of the gravitino 
mass^. For example for 777,3/2 of the order of 1-1000 GeV one obtains ~ 0.1 

for Tr ~ 10^ - 10^^ GeV, with ~ 1 TeV. Even with a high value of Tr there is 
no gravitino problem, since the next-to-LSP decays to standard model particles much 
earlier than BBN epoch via R-parity breaking interactions. 

Let us now show that If/^i^P ~ 10~^® — 10~^^ in the /^.t^SSM. We can easily make an 
estimation. For a 2 x 2 matrix, 

a c 



c b 



(14) 



^Let us recall that there is a lower limit of 1.2 keV on the mass of (warm) dark matter particles 
from Lyman a forest [26]. 
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the mixing angle is given by tan 2^ = 2c/ {a — b). In our case (see Eq. (5)) c ~ 
QiU ~ 10^^ GeV (represents the mixing of Bino and left handed neutrino), a ~ 1 TeV 
(represents the Bino mass Mi), and 6 = 0. Thus one obtains tan 26' ~ 10~^, implying 
sin^ ~ ^ ~ 10-^. This gives jK^^p ~ lO^^^ More general, ^ ~ lO'^ - lO'^, 

giving rise to 

io-''< if^7-r<io-''. (15) 

In order to confirm this estimation we have performed a scan of the low-energy 
parameter space of the model discussed in Sect. 2, over the following ranges: 
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Concerning the rest of the soft parameters, we will take for simplicity in the compu- 
tation rriQ -c jc gc = 1 TeV, Ay_^d,e — 1 TeV, Ax = —A^ = — 2^^ = 1 TeV, and for 
the other gaugino masses we will use the GUT relations. Although this is not a full 
exploration of the parameter space, which is beyond the scope of this work, it gives a 
fair estimation of the representative values for If/^yi^p. The results are shown in Fig. 1. 
The black points there correspond to regions of the parameter space where the current 
data on neutrino masses and mixing angles are reproduced (where we are using the 
allowed 3cr ranges discussed in [27]). In addition, these regions avoid false minima and 
tachyons, as well as fulfil the Landau pole constraint, following the lines discussed in 
[2, 11]. Typically, the mass of the lightest neutralino is above 20 GeV, and since the 
gravitino mass in this model is constrained to be below that value, as we will see in 
the next section, the gravitino can be used as the LSP. 

In principle, we could conclude that the range 10"^^ < |f^7i^P ^ 5 x 10"^'^ is specially 
favoured. Nevertheless, despite that we see only a few solutions for jf/^iyp < lO""*^^, 
looking at Eq. (11) we could infer that values close to zero would be achievable through 
a cancellation of the Bino and Wino contribution. Therefore we consider that a good 
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Figure 1: It/^z/p versus the lightest neutrino mass. 

estimation for the lower bound of \U^u\^ without much fine tunning is 10~^^. On the 
other hand, one could get values \U^v\'^ > 10~^^ in the regime of degenerated neutrinos, 
where larger values of the lightest neutrino mass than the ones shown in Fig. 1 are 
required (i.e. allowing an exploration for larger values of Y^- and the VEVs z/j). In this 
sense, one could use the (conservative) range written in Eq. (15). This is what we will 
do in the next section. 

4 Gamma rays from gravitino decay 

Since in R-parity breaking models the gravitino decays producing a monochromatic 
photon with an energy 7713/2/2, one can try to extract constraints on the parameter 
space from gamma-ray observations [14,28]. Actually, model independent constraints 
on late dark matter decays using the gamma rays were studied in [29] (see also also [30] 
for the case of neutrino production). There, the decaying dark matter was constrained 
using the gamma-ray line emission limits from the galactic center region obtained 
with the SPI spectrometer on INTEGRAL satellite, and the isotropic diffuse photon 
background as determined from SPI, COMPTEL and EGRET data. These constraints 
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Figure 2: Constraints on lifetime versus mass for a decaying dark matter particle. The 
region below the magenta solid line is excluded by gamma-ray observations [29]. The 
region below the green dashed (blue dotted) line will be checked by FERMI. Black solid 
lines correspond to the predictions of the /iz/SSM for several representatives values of 

|?7;^,|2 = 10-16 _ 10-12. 

are shown with a magenta line in Fig. 2, where a conservative non-singular profile at 
the galactic center is used. 

On the other hand, the FERMI satellite [23] launched in June 2008 is able to 
measure gamma rays with energies between 0.1 and 300 GeV. We also show in Fig. 2 
the detectability of FERMI in the 'annulus' and 'high latitude' regions following the 
work in [18]. Below the lines, FERMI will be able to detect the signal from decaying 
dark matter. Obviously, no signal means that the region would be excluded and FERMI 
would have been used to constrain the decay of dark matter [18]. 

Finally, we show in the figure with black solid lines the values of the parameters 
predicted by the /xz/SSM using Eq. (12), for several representative values of If/^^^p. 
We can see that values of the gravitino mass larger than 20 GeV are disfavored in 
this model by the isotropic diffuse photon background observations (magenta line). In 
addition, FERMI will be able to check important regions of the parameter space with 
gravitino mass between 0.1 — 20 GeV and jf/^j/p = lO^^^ — lO^^^ (those below the green 
line) . 
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Let us now discuss in more detail what kind of signal is expected to be observed by 
FERMI if the gravitino lifetime and mass in the /xi^SSM (black solid lines) correspond 
to a point below the green line in Fig. 2 

As it is well known, there are two sources for a diffuse background from dark matter 
decay. One is the cosmological diffuse gamma ray coming from extragalactic regions, 
and the other is the one coming from the halo of our galaxy. 

The photons from cosmological distances are red-shifted during their journey to the 
observer and the isotropic extragalactic flux turns out to be [14, 18, 28] 

= A^,^ { 1 ^.{^)X' f — )" e (i - ^) . m 

dE moM \^ V"^DM/ J \mBMj \ mBuJ 

with 

- , = 2.11 X io-( s str)-' {^y i^y. m 

Here k = Qa/^m ^ 3 with Qa + Qm = 1, Pc = 1-05 h'^ x IQ-^GeVcm-^, Ho = 
h 100 km s^^ Mpc^^ with h = 0.73 [31], and Tdm and ttt-dm are the lifetime and 
mass of the dark matter particle, respectively. We take the dark matter density as 
^DMh'^ = 0.1. 

On the other hand, the photon flux from the galactic halo shows an anisotropic 
sharp line. For decaying dark matter this is given by 

«^ = A,^^J-6 (l - — ) . (19) 
dE rriBM \ rnoM J 

with 



Ahalo = / Phalo{l)dl , (20) 

47rTDM"T.DM Jlos 

where the halo dark matter density is integrated along the line of sight, and we will 
use a NFW profile [32] 

^""^^'^ ^ r/r/+r/r.y ' (^1) 
where we take ph — 0.33 GeV/ cm^ = 0.6 x lO^Pc, fc — 20 kpc, and r is the distance 
from the center of the galaxy. The latter can be re-expressed using the distance from 
the Sun, s, in units of Rq — 8.5 kpc (the distance between the Sun and the galactic 
center) and the galactic coordinates, the longitude, I, and the latitude, b, as 

y{s, b, I) = Ry[{s - cos b cos if + (1 - cos^ b cos^ /)] . (22) 
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E [GeV] E [GeV] 

(a) (b) 

Figure 3: Expected gamma-ray spectrum for an example of gravitino dark matter 
decay in the mid-latitude range (10° < |6| < 20°) in the /xt^SSM with m3/2 = 3.5 GeV 
and (a) \U^^f = 8.8 x lO'^^ corresponding to ry2 = 10^^ s, (b) {U^y^f = 1.7 x 10"^^ 
corresponding to r3/2 = 5 x 10^^ s. The green dashed, magenta solid, and black solid 
lines correspond to the diffuse extragalactic gamma ray flux, the gamma-ray flux from 
the halo, and to the conventional background, respectively. The total gamma-ray flux 
is shown with red solid lines. The blue solid lines are explained in the note added in 
Sect. 6. 
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Figure 4: The same as in Fig. 3 but for 7713/2 = 10 GeV and (a) \U^uf = 3.8 x 10~ 
corresponding to r3/2 = 10^'' s, (b) \U^uf = 7.6 x 10^"'^'^ corresponding to r3/2 = 5 x 10' 
s. 



It is worth noticing here that when computing above the gamma-ray fluxes, the 
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effects of attenuation of the flux in the interstellar or the intergalactic medium have 
been neglected. In our mass range of the gravitino, the flux from the decay of the 
gravitino dark matter is made of photons and neutrinos, thus we might expect the 
attenuation of the gamma-ray flux by pair production. Nevertheless, for our case with 
less than 10 GeV gamma-ray flux the attenuation is suppressed both in the galactic 
and extragalactic medium [33], and can therefore be safely neglected. 

Let us now compute with the above formula, as an example, the expected diffuse 
gamma-ray emission in the mid-latitude range (10° < |6| < 20°), which is being 
analized by FERMI [34], for the case of gravitino dark matter. Let us assume for 
instance a value of m3/2 = 3.5 GeV and \U^u^ = 8.8 x 10^^^(1.7 x 10"^^) in the 
/iz/SSM, corresponding to r3/2 = 10^'^ (5 x 10^^) s, using Fig. 2. We convolve the 
signal with a Gaussian distribution with the energy resolution AE/E = 0.08, between 
= 1 — 10 GeV, following [23], and then we average the halo signal over the region for 
the mid-latitude range mentioned above. 

The results for the two examples are shown in Fig. 3. There, the green dashed 
line corresponds to the diffuse extragalactic gamma ray ffux, and the magenta solid 
line corresponds to the gamma-ray flux from the halo. The black solid lines represent 
the background including the diffuse galactic emission model from GALPROP [35], 
and point source and isotropic contributions [34] . The systematic uncertainties for the 
latter generate the band shown within the two black lines. 

The total gamma-ray ffux, including background, extragalactic, and line signal, is 
shown with red solid lines. We can see that the sharp line signal associated to an 
energy half of the gravitino mass, dominates the extragalactic signal and can be a 
direct measurement (or exclusion) in the FERMI gamma ray observation. We could 
also use the EGRET data [36] to constrain the parameter space of the model. Although 
the data beyond 1 GeV are controversial (even instrumental effects might be a possible 
explanation for the observed excess [37]), the line obtained in the example of Fig. 3a 
is very sharp and could be discarded when compared with the spectrum of EGRET. 

As another example, we show in Fig. 4 the case of 7713/2 = 10 GeV for the same 
values of lifetimes as in Fig. 3. 
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Figure 5: Constraints on lifetime versus mass for gravitino dark matter in the /ii^SSM. 
The region below the magenta solid line is excluded by several gamma-ray observations 
[29]. The region below the red solid line is disfavoured by FERMI. Black solid lines 
correspond to the predictions of the /xi^SSM for several representatives values of lU^ul"^ = 

10-16 _ iQ-12_ 

5 Conclusions 

We have discussed the possibility of gravitino dark matter in the /iz/SSM, where R- 
parity is broken and therefore the LSP is unstable. For the gravitino mass of the 
order of GeV, the lifetime of the gravitino is much longer than the age of the Universe 
and the observed relic density can be explained well by the thermal production of 
gravitinos after inflation. As a consequence, the gravitino can be a good candidate for 
dark matter. 

We have also studied the prospects for detecting gamma rays from decaying grav- 
itinos. If the gravitino explains the whole dark matter component, the gravitino mass 
larger than 20 GeV is disfavored by the isotropic diffuse photon background measure- 
ments. Nevertheless, a gravitino with a mass range between 0.1 — 20 GeV gives rise to 
a signal that might be observed by the FERMI satellite. In this way important regions 
of the parameter space of the /iz/SSM can be checked. 

6 Note added 

After completion of the current work, the FERMI experiment reported 5-month mea- 
surements of the diffuse gamma-ray emission in the mid-latitude range [38]. We have 
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added in Figs. 3 and 4 blue solid lines corresponding to these Fermi LAT data with 
systematic uncertainties. These turn out to be consistent with the background model, 
implying that the sharp lines obtained in the examples of Figs. 3a and 4a have not 
been observed. Taking these results into account, we have summarized in Fig. 5 the 
constraints on lifetime versus mass for the /luSSM. Values of the gravitino mass larger 
than 10 GeV are now disfavored, as well as lifetimes smaller than about 3 to 5 x 10^^ s. 
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